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Abstract 


The  first  algorithm  for  solving  generalized  linear  programs  was  given  by 
— George  Dantzig.  His  algorithm  assumes  that  a  basic  feasible  solution  of  the 

generalized  linear  program  to  be  solved  exists  and  is  given.  If  the  initial  basic  fea¬ 
sible  solution  is  non-degenerate,  then  his  algorithm  is  guaranteed  to  converge.  The 
purpose  of  this  paper  is  to  show  how  to  find  an  initial  basic  feasible  (possibly  de¬ 
generate)  solution  of  a  generalized  linear  program  by  applying  the  same  algorithm 
to  a^Rphase-one”*problem  without  requiring  that  the  initial  basic  feasible  solution 
to  the  latter  be  non-degenerate.  ^ 

Key  Words:  linear  programming,  generalized  linear  programming,  duality. 
Abbreviated  Title:  The  feasibility  of  a  GLP 


Generalized  linear  programming  is  a  natural  extension  of  linear  programming. 
The  cost  coefficients  and  the  columns  of  the  constraint  matrix  of  a  generalized  linear 
program  may  be  freely  chosen  from  a  compact  convex  set.  Specifically,  a  generalized 
linear  program  is  defined  as 

(GLP)  minimize  cy 
subject  to 

py  =  q 

( c,pT)T  may  be  freely  chosen  from  T 1  x  Tm 
y  >0 

where  y  >  0  is  a  real  number,  q  is  in  Rm ,  T1  is  a  closed  interval  in  R1 ,  and  Tm  is 
a  compact  convex  set  in  Rm.  - — -  _  '  . 
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Dantzig’s  algorithm  for  solving  generalized  linear  programs  starts  with  a  basic 
feasible  solution  of  the  problem  to  be  solved,  i.e.,  assumes  that  one  has  at  hand  m 
independent  vectors  p1 ,  p2,  . . . ,  pm  such  that  p*  £  Tm  and  values  of  scalars  y*  >  0 

such  that  YZLiP'Vi  =9*s  non-negative.  At  iteration  k,  a  new  vector  (c*,p*  )T  £ 
T 1  x  Tm  is  generated  with  the  most  negative  reduced  cost  ck  —  irpk  <  c  —  irp  for 
all  (c,pT)T  £  T1  x  Tm,  where  Trp*  =  c*  for  basic  pJi.  Then  p*  is  added  to  the 
constraint  matrix  of  the  k- th  restricted  master  program  with  cost  coefficient  ck .  If 
the  initial  basic  feasible  solution  is  non- degenerate,  then  the  optimal  solutions  of  all 
restricted  master  programs  form  a  sequence  of  feasible  solutions  of  (GLP)  on  which 
the  objective  function  tends  to  the  value  of  the  program  (GLP)  (see,  e.g.,  Dantzig 
(1963),  chapters  22  and  24). 

The  purpose  of  this  paper  is  to  show  how  to  find  an  initial  basic  feasible  (pos¬ 
sibly  degenerate)  solution  of  (GLP)  by  Dantzig’s  algorithm.  Here  the  idea  is  same 
as  phase-one  of  the  simplex  method.  We  introduce  artificial  variables  and  minimize 
the  summation  of  all  artificial  variables.  That  is,  we  solve  the  phase-one  problem: 

(LP)  minimize  x, ,• 

subject  to 

Py  +  Y,?=i  e'xi  =  Q 
p  may  be  freely  chosen  from  Tm 

y  >  0,  X{  >  0  for  i  =  1,  . . . ,  m, 

where  e*  is  the  z'-th  unit  vector  in  Rm.  We  assume,  without  loss  of  generality,  that 
qi  >  0  for  i  —  1,  ...  m,  and  g  /  0.  In  addition,  we  assume  that  the  origin  is 
not  contained  in  Tm.  Let  u(LP)  denote  the  optimal  value  of  the  objective  function 
of  program  (LP).  Superscripts  on  vectors  denote  different  vectors,  while  subscripts 
on  vectors  denote  different  components.  To  avoid  confusion,  we  state  the  steps  for 
finding  an  initial  basic  feasible  solution  of  (GLP)  explicitely. 
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Initialize. 

Let  k  :=  0; 

let  the  first  restricted  master  linear  program,  LP(0),  be 
minimize  x, 

subject  to 

,=i  e*Xj  =  q 

Xi  >  0  for  i  =  1,  . . . ,  m. 

(Note  that  v(LP(0))  =  >  0  since  g  >  0  and  q  ^  0.) 

Step  1. 

If  u(LP(k))  =  0,  then  a  basic  feasible  solution  of  (GLP)  is  found,  stop. 

Else,  go  to  Step  2. 

Step  2. 

Let  Trk  be  an  optimal  dual  solution  of  LP(k); 
find  a  pk+l  such  that  nkpk+i  =  max{7rfcp  :  p  €  Tm); 
if  7 rfcpfc+1  <  0,  then  (GLP)  is  infeasible,  stop. 

Else,  go  to  Step  3. 

Step  3. 

Let  LP(k+l)  be  the  linear  program 
minimize  2^1=1  xi 
subject  to 

Ei=i  P'Vi  +  £.=1  e’*»  =  9 

y,  >  0,  *  =  1,  . . . ,  k  +  1,  Xi  >  0,  i  =  1,  . . . ,  nv, 

solve  LP(k+l); 
k  :=  k  +  1; 
go  to  Step  1. 

Theorem.  Suppose  that  the  origin  is  not  contained  in  Tm.  Then  (GLP)  is 
feasible  if  and  only  if  t/(LP(k))  =  0  for  some  k  or  v(LP(k))  >  0  for  all  k  and 
limt—oo  i>(LP(k))  =  0. 
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Proof.  Suppose  first  that  t>(LP(k))  >  0  for  all  k  and  lim*- .00  v(LP(k))  =  0.  Let 
(l/i  1  •  •  •  1  Mkix  1  > •  •  •  ■> xm)  be  an  optimal  solution  of  LP(k).  Let  A*  be  the  index  set 
ordered  from  low  indicies  i  to  high  indecies  such  that  yk  >  0  if  and  only  if  i  G  A* .  If 
we  solve  LP(k)  by  the  simplex  method,  then  A*  contains  no  more  than  m  elements. 
Let  A*  =  {**, A|tl}  where  |At|  denotes  the  cardinal  number  of  A *,  and 
let  Yk  =  (yk.  , ...  ,y*  ,0,...  ,0)  G  Rm,  that  is,  the  first  |At|  components  of 

*1  *1^*1 

Yk  are  yt*  , ...  , ykk  and  the  rest  are  zeros  in  the  case  |A*|  <  m.  We  claim 

that  {Yk  G  Rm  :  k  =  1,2,...}  has  a  cluster  point.  Indeed  let  \k  = 

for  all  k  =  1,  2, _ Since  A*  =  0  implies  t>(LP(k))  =  ]>^i  q%  >  0,  we  know  that 

A*  =  0  can  only  happen  for  a  finite  number  of  k  because  an  infinite  number  of  times 
would  imply  limJ_00  u(LP(kj))  =  YaLiQ*  >  0  for  a  subsequence,  contrary  to  the 
hypothesis  lim*— v(LP(k))  =  0.  Therefore,  A*  >  0  for  all  k  >  ko.  Now  let  sk  = 
]Cf=i  p'yf/A*  G  Tm  and  uk  =  e'xk  for  all  k  >  k0.  As  lim*— t>(LP(k))  = 
0,  we  have  lim*— oo  uk  =  0  and  lim*— oo  Xksk  —  q.  Since  Tm  is  compact,  there 
exists  a  subsequence  sk>  such  that  limj_0Osfc>  =  s*.  We  know  that  s*  ^  0  since, 
by  hypothesis,  the  compact  set  Tm  does  not  contain  the  origin.  It  follows  that 
limJ_oo^*>  =  ||<?||  •  ||s*||-1  and  thus  \k>  is  bounded.  Therefore,  the  sequence  Yk 
has  a  cluster  point.  Let  Ak  =  (p'*i , . . .  ,p**iA*i , 0, . . . ,  0)  G  RmXm  for  all  k  >  0.  As 
Tm  is  compact,  for  a  subsequence  kj  we  have 

lim  Ak>  =  A*,  lim  Yk>  =  Y\  and  lira  {Ak*Yki  +  uk* )  =  A*Y *  =  q. 

j—+oo  j— ►  OO  j— *  00 


Now  let  p1*,  . . . ,  pr*  be  the  nonzero  columns  of  A*,  yi*,  . . . ,  yOT*  be  the  components 
of  Y*,  and  choose  p(r+1^*,  . . . ,  pu  G  Tm  such  that  ran^p1*, . . .  ,p‘*)  =  m.  Then 
Vi  =  Him  for  i  —  1,  . . . ,  r,  and  y<  =  0  for  i  =  r  +  1,  . . . ,  t  is  a  feasible  solution  of  the 
system: 

plmyi  +---  +  Pt*yt  =  q 

y»  >  0,i  =  ' 


and  thus  we  can  find  a  basic  feasible  solution  of  (1)  by  the  simplex  method.  Of 
course,  this  basic  feasible  solution  is  also  a  basic  feasible  solution  for  (GLP). 

It  remains  to  show  that  lim*- 00  v(LP(k))  >  0  implies  (GLP)  infeasible.  Since 
Trk  is  an  optimal  dual  solution  of  LP(k), 


7rk  <  1  for  i  =  1,. . .  ,m, 


(2) 


and 

m 

t kq  =  xk  >  0  for  all  k  =  1, 2, . . .  .  (3) 

1=1 
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Since  q  >  0,  (2)  and  (3)  tell  us  that  the  sequence  nk  is  bounded.  The  boundedness 
of  7r*  implies  lim^oo  v(LP(k))  =  u(LP)  (see,  e.g.,  Dantzig  (1963),  chapter  24).  It 
follows  that  u(LP)  >  0  and  (GLP)  is  infeasible.  | 

Remarks.  Given  any  (GLP)  with  0  ^  q  >  0,  one  can  always  generate  a  sequence 
of  linear  programs  LP(k).  Since  {7rfc  :  k  =s  1,2, . . .  }  is  bounded,  limfc_00  u(LP(k))  = 
v(LP).  Therefore,  lim*..,,*,  v(LP(k))  >  0  always  implies  that  (GLP)  is  infeasible. 
While  in  the  case  limfc_00  v(LP(k))  =  0,  if  limsupt^^  ||.sfc||  >  0,  then  v(LP)  =  0 
can  be  attained,  which  implies  that  (GLP)  is  feasible. 
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